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ABSTRACT 

In  the  present  study,  nine  empirical  methods  for  calculating  daily  reference  evapo  transpiration  (ET0)  namely, 
Blaney-Criddle,  Jensen-Haise  and  Hargreaves  (temperature  based),  Priestley-Taylor,  Radiation  and  Makkink 
(radiation  based),  Modified  Penman  (physically  based),  Pan  Evaporation  and  Christiansen  (pan  evaporation  based) 
methods  have  been  evaluated  with  respect  to  FAO-56  Penman-Monteith(PM)  method  for  estimating  daily  ET0  in  the 
semi-arid  Nellore  region  of  Andhra  Pradesh,  India.  Data  was  collected  from  the  India  Meteorological  Department  (IMD), 
Pune.  The  evaluation  is  based  on  performance  criteria  namely,  Root  Mean  Square  Error  (RMSE),  Coefficient  of 
Determination  (R2)  and  Efficiency  Coefficient  (EC).  The  relationships  between  PM  method  and  the  other  methods  were 
developed  to  obtain  daily  ET0  estimates  comparable  with  PM  method.  The  ET0  equations  were  then  recalibrated  with 
respect  to  PM  method  for  improving  their  daily  ET0  estimation  capability  in  the  region  selected  for  the  present  study. 
The  recalibrated  Modified  Penman  and  Blaney-Criddle  methods  showed  satisfactory  performance  in  the  daily  ET0 
estimation.  However,  the  recalibrated  Blaney-Criddle  method  may  be  adopted  because  of  its  simpler  data  requirements 
with  reasonable  degree  of  accuracy. 

KEYWORDS:  Reference  Evapotranspiration,  Recalibration,  Performance  Evaluation 
INTRODUCTION 

A  reliable  estimation  of  Evapotranspiration  (ET)  is  of  critical  importance  in  irrigation  system  design,  crop  yield 
simulation  and  water  resources  planning  and  management.  Field  measurement  of  evapotranspiration  is  rarely  available  and 
actual  crop  evapotranspiration  (ETC)  is  usually  calculated  from  reference  evapotranspiration  (ET0)  using  the  crop  factor 
method,  which  consists  of  multiplying  ET0  with  crop  coefficients  (Kc)  to  obtain  ETC  (i.e.,  ETC  =  ET0  x  Kc).  Several  reports 
on  the  estimation  of  Kc  are  available. Allen  et  al.  (1998)[2]  and  Jensen  et  al.  (1990)'51have  reported  crop  coefficients  for 
many  crops.  These  values  are  commonly  used  in  places  where  the  local  data  is  not  available. 

It  is  desirable  to  have  a  method  that  estimates  reasonably  the  reference  Evapotranspiration  (ET0).  According  to  the 
Food  and  Agricultural  Organization  (FAO),  FAO-56  Penman-Monteith  (PM)  method,  that  requires  numerous  climatic 
parameters,  achieves  better  agreement  with  the  lysimeter  ET0  measurements  compared  to  all  other  known  methods. 
However,  under  limited  climatic  data  availability  conditions,  the  simple  empirical  methods  yielding  results  comparable 
with  PM  ET0  may  be  selected  at  regional  level  for  reasonable  estimation  of  ET0. 

Irmak  et  al.(2003)[41,  Berengena  and  Gavilan(2005)[3],  Alkaeed  et  al.(2006)[l],  Suleiman  and  Hoogenboom 
(2007)[S1,  Trajkovic  and  Stojnic  (2008)[9]  have  evaluated,  compared  and  tested  the  applicability  of  various  ET0  equations 
for  different  regions. 
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The  present  study  reports  the  performance  evaluation  of  commonly  used  nine  ET0  estimation  methods  based  on 
their  accuracy  of  estimation  and  these  methods  are  recalibrated  with  PM  method  for  Nellore  region  of  Andhra  Pradesh 

MATERIALS  AND  METHODS 

Nellore  region,  located  in  Nellore  district  of  Andhra  Pradesh,  India,  with  global  coordinates  of  14°  22  N  latitude 
and  79°  59  E  longitudes,  has  been  chosen  as  the  study  area.  Daily  meteorological  data  at  the  region  for  the  period 
1983-2003  was  collected  from  India  Meteorological  Department  (IMD),  Pune.  A  part  of  the  data  (1983-1997)  was  used  for 
developing  recalibrated  equations,  while  the  rest  of  the  data  (1998-2003)  was  used  to  verify  the  performance  of  the 
recalibrated  equations.  The  details  of  the  methods  selected  for  the  present  study  are  presented  in  Table  1. 


Table  1:  Details  of  Reference  Evapotranspiration  Estimation  Methods 
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PERFORMANCE  EVALUATION  CRITERIA 


The  performance  evaluation  criteria  used  in  the  present  study  are  the  coefficient  of  determination  (R~),  the  root 
mean  square  error  (RMSE),  systematic  RMSE,  unsystematic  RMSE  and  the  efficiency  coefficient  (EC). 

Coefficient  of  Determination  (R2) 

It  is  the  square  of  the  correlation  coefficient  (R)  and  the  correlation  coefficient  is  expressed  as 


^  (Oi-oXPi-p) 


R  = 


-1 1/2 


Where  O  and  P  are  observed  and  estimated  values,  ^  and  P  are  the  means  of  observed  and  estimated  values  and 
n  is  the  number  of  observations.  It  measures  the  degree  of  association  between  the  observed  and  estimated  values  and 
indicates  the  relative  assessment  of  the  model  performance  in  dimensionless  measure. 

Root  Mean  Square  Error  (RMSE) 

It  yields  the  residual  error  in  terms  of  the  mean  square  error  and  is  expressed  as  (Yu  et  al.,  1994)[101 


RMSE  =  "l,J=L 


n 

Systematic  RMSE  (RMSES) 

It  measures  the  room  available  for  local  adjustment.  It  is  expressed  as 


RMSE 


1 


Z(A-*,)2 


Where  P'    a  +     ,  a  and  b  are  the  liner  regression  coefficients 

Unsystematic  RMSE  (RMSEU) 

It  shows  the  noise  level  in  the  model  and  is  a  measure  of  scatter  about  the  regression  line  and  potential  accuracy. 
It  is  expressed  as 


RMSE,  = 


Z  (Pi- fry 


1=1 
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Efficiency  Coefficient  (EC) 

It  is  used  to  assess  the  performance  of  different  models  (Nash  and  Sutcliffe,  1970)m.  It  is  a  better  choice  than 
RMSE  statistic  when  the  calibration  and  verification  periods  have  different  lengths  (Liang  et  al.,  1994)[61.  It  measures 
directly  the  ability  of  the  model  to  reproduce  the  observed  values  and  is  expressed  as 

±(o,-p,Y 

EC  =  1  -   

A  value  OF  EC  of  90%  generally  indicates  a  very  satisfactory  model  performance  while  a  value  in  the  range  of 
80-90%  indicates  a  fairly  good  model.  Values  of  EC  in  the  range  60-80%  would  indicate  an  unsatisfactory  model  fit. 

RESULTS  AND  DISCUSSIONS 

ET0  was  estimated  using  the  climatic  data  for  different  methods  with  original  values  of  empirical  coefficients. 
The  mean  daily  values  were  compared  with  those  estimated  by  PM  method.  The  percentage  deviations  with  reference  to 
PM  method  are  shown  in  Table  2.  The  positive  deviation  represents  overestimation  and  negative  deviation  indicates 
underestimation  of  ET0  values.  It  may  be  observed  that  the  percentage  deviation  between  PM  method  and  other  methods 
varied  between  -14.9%  and  55.3%.  RA  method  estimated  ET0  with  largest  deviation  followed  by  JH  method  and,  PT,  CS, 
HR  and  BC  methods  with  minimum  deviation.  The  performance  indicators  of  the  methods  with  original  coefficients  are 
presented  in  Table  3.  The  relatively  more  unsystematic  RMSE  components  with  the  ET0  estimation  methods  except  MP 
and  BC  methods  indicate  more  noise  level  in  the  methods  and  scatter  about  the  regression  line. 

The  temperature,  radiation,  physically  and  pan  evaporation  based  methods  selected  for  the  present  study  were 
recalibrated  with  respect  to  PM  method  as  presented  in  Table  4.  The  performance  indicators  of  these  empirical  models  with 
original  and  recalibrated  coefficients  in  the  estimation  of  ET0  are  given  in  Table  5.  The  improved  R2,  EC  and  reduced 
RMSE  (Table  5)  indicate  the  closeness  of  estimated  daily  ET0  values  and  thereby  reflect  the  appropriateness  of 
recalibration.  An  improvement  in  the  performance  of  ET0  estimation  methods  with  recalibrated  coefficients  over  these 
methods  with  original  coefficients,  in  general,  has  been  observed  (Table  5).  A  significant  improvement  has  been  found  in 
case  of  recalibrated  MP  and  BC  methods.  However,  out  of  these  methods,  recalibrated  BC  method  may  be  adopted  in  the 
reasonable  daily  ET0  estimation  in  the  region  because  of  simpler  data  requirements.  The  scatter  plots  as  shown  in  Figure  1 
&  2  also  depict  similar  observations. 


Table  2:  Percentage  Deviations  in  the  Estimated  Mean  Daily  Reference 
Evapotranspiration  with  Original  Coefficients 
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Table  3:  Performance  Indicators  of  Various  Methods  with  Original  Coefficients  against  PMM 


Method 

Slope(m) 

Intercept(c) 

R2 

RMSE 
(mm) 

RMSE 
s  (mm) 

RMSEU 
(mm) 

EC  (%) 

BC 
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JH 
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HR 
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0.90 
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0.69 
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PT 
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0.2826 
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0.86 

0.52 

0.68 
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RA 

0.5547 

0.6238 
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0.92 

0.60 

0.70 
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MK 

1.0149 
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0.5740 

0.92 

0.60 

0.70 
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MP 
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0.26 

0.05 

0.26 
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PE 

0.5523 

2.3673 

0.3780 

1.11 

0.88 

0.68 
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CS 

0.8921 

0.3792 
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0.50 

0.18 

0.47 
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Table  4:  ET„  Estimation  Methods  with  Original  and  Recalibrated  Coefficients 
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Table  5:  Performance  Evaluation  of  ET0  Estimation  Methods  with  Original 
and  Recalibrated  Coefficients  against  PM  Method 
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Figure  1:  Scatter  Plots  of  Daily  ET0  Estimated  by  Various  Methods  with  Original 
Coefficients  against  ET0  Estimated  Using  PM  Method 
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Figure  2:  Scatter  Plots  of  Daily  ET0  Estimated  by  Various  Methods  with  Recalibrated  Coefficients  against  ET0 

Estimated  using  PM  Method  during  Testing  Period 
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CONCLUSIONS 

The  Blaney-Criddle,  Jensen-Haise  and  Hargreaves  (temperature  based),  Priestley-Taylor,  Radiation  and  Makkink 
(radiation  based),  Modified  Penman  (physically  based),  Pan  Evaporation  and  Christiansen  (pan  evaporation  based) 
reference  evapotranspiration  (ET0)  methods  have  been  recalibrated  with  respect  toFAO-56  Penman-Monteith  method  and 
their  performance  in  the  daily  ET0estimation  in  the  Nellore  region  of  Andhra  Pradesh  has  been  evaluated.  All  these  ET0 
estimation  methods,  in  general,  showed  an  improved  performance  with  recalibrated  coefficients.  But,  the  recalibrated 
Modified  Penman  method  and  recalibrated  Blaney-Criddle  method  performed  well  in  the  daily  ET0  estimation.  However, 
recalibrated  Blaney-Criddle  method  may  be  applied  for  the  reasonable  estimation  of  daily  ET0  in  the  region  because  of 
simpler  data  requirements. 
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